Orthogonality constraints and entropy in the SO(5)-Theory of 

HighT c -Superconductivity 

Franz J. Wegner 
Institut fur Theoretische Physik 
Ruprecht-Karls-Universitat Heidelberg, Germany 

February 1, 2008 

Abstract 

S.C. Zhang has put forward the idea that high-temperature-superconductors can be described 
in the framework of an SO(5)-symmetric theory in which the three components of the antifer- 
romagnetic order-parameter and the two components of the two-particle condensate form a five- 
component order-parameter with SO(5) symmetry. Interactions small in comparison to this strong 
interaction introduce anisotropies into the SO(5)-space and determine whether it is favorable for 
the system to be superconducting or antiferromagnetic. Here the view is expressed that Zhang's 
derivation of the effective interaction V a a based on his Hamiltonian _ff a is not correct. However, 
the orthogonality constraints introduced several pages after this 'derivation' give the key to an 
effective interaction very similar to that given by Zhang. It is shown that the orthogonality con- 
straints are not rigorous constraints, but they maximize the entropy at finite temperature. If the 
interaction drives the ground-state to the largest possible eigenvalues of the operators under con- 
sideration (antiferromagnetic ordering, superconducting condensate, etc.), then the orthogonality 
constraints are obeyed by the ground-state, too. 

Key words: Superconductivity; Strongly correlated electrons 

1 Introduction 

In a recent paper |jj Zhang has proposed the idea that the components of the electron-pair conden- 
sate and those of the antiferromagnetic order-parameter form a five-component order parameter in 
an approximately SO(5)-invariant theory. Anisotropies small in comparison to the SO(5)-invariant 
interaction break this symmetry. They depend in particular on the chemical potential and thus on 
doping. This allows him to desribe the transition from antiferromagnetism in the half-filled system to 
superconductivity under moderately weak doping. 

This phenomcnological picture has been tested by numerical calculations for small Hubbard- 
systems ^ [| and t- J- models [|| |), for which good agreement with the predictions from Zhang's 
theory are obtained. It has also been shown that spin-ladders observe SO(5)-symmetry j|, || or even 
SO(8)-symmetry 0. 

Returning to higher-dimensional systems, it is pointed out here that Zhang's derivation of the effec- 
tive interaction V c s based on his Hamiltonian H a seems to be not correct. (Actually we struggled for 
quite a while trying to understand his derivation). However, the orthogonality constraint introduced 
several pages after this 'derivation' gives the key to an effective interaction very similar to that given 
by Zhang. After having established this in sect. 2 it is shown in sect. 3, that the entropy contains 
squares of these constraint terms. Thus the orthogonality constraints are not strict requirements, but 
their fulfillment maximizes the entropy and lowers the free energy. Finally in sect. 4 it is argued, that 
also in the ground-state these constraints are likely to hold. 
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2 Effective Interaction and the Orthogonality Constraint 



Zhang introduces a five-component order-parameter n i; i = 1..5. The operators associated with these 
components are given explicitcly in the appendix. Two components, n\ and n 5 , describe the real part 
and the imaginary part of the superconducting condensate, the three other components ri2, n^, and 
n4, are the three Cartesian components of the staggered magnetization. He assumes that a strong 
interaction which may be described by a Ginzburg-Landau interaction, leads to a symmetry breaking 
below a critical temperature T c and has an SO(5) symmetry, thus allowing both superconducting and 
antiferromagnetic order. This leading interaction is a function of J2 i nf only and does not give any 
preference to superconducting or antiferromagnetic order. 

In addition there is a weaker interaction H a which introduces anisotropy into the system 



u ^1,5 . ^1.2 + ^1.3 + ^1.4 + ^2.5 + ^3.5 + ^4.5 . -^2,3 + ^2.4 + ^3,4 9 , 2, 2 2\ or 

B.= We + ^ + ^ -(n 2+ n 3 + n 4 )-2^ 1 ,, 

(1) 

The L a j, are operators bilinear in the electron creation and annihilation operators. They are given 
cxplicitely in appendix A. [i is the chemical potential. The operators L obey L b ^ a = — L a ^ and the 
commutator relations 

[L a ,b, L c .d] = -i5b, c L a ,d + iS a ,cLb,d + ih,dL a>c - i5 a>d L b ^ c (2) 

of the orthogonal group SO(5), moreover the vector n is rotated by L, 

[L a>b , n c ] = -i5 b , c n a + iS a ^ c n c . (3) 

The derivative terms of Zhang's interaction are left out here, since only the global state is investigated. 
The effective potential is determined as the minimum of the Hamiltonian for given order-parameter 
m, ... ri5, which Zhang normalizes to 

E«? = 1 - (4) 

i 

Since one considers the interaction for a macroscopic system, all the components of L and n can be 
considered as classical quantities and the obvious minimum is obtained for L\*, = 2^x c , whereas all 
other components of L vanish, which yields what I call V^g lv 

V c t v (n) = -2^ Xc - 9 -{n 2 2 + nj+ n|). (5) 

It would mean that the only effect of the chemical potential is a lowering of the energy of the system, 
but it would have no effect on an anisotropy in the order-parameter space. This obviously is at 
variance with Zhang's claim for the effective interaction 

K Z ff hanS H = -2M 2 («i + n§)[Xc(n? + ng) + X M + n\ + n\)\ - \{n\ + n§ + n\). (6) 

An indication that Zhang's Lagrangian 

^ = Y,^f^ + V{n) (7) 

a<b 

^ab = n a (d T n h - iB bc n c ) - (a -> b) (8) 

is not equivalent to iJ a , is obvious from the fact that in Zhang's Lagrangian the only allowed rotations 
are in the plain spanned by n and d T n — iBn, whereas it does not permit a rotation around a 
second perpendicular plain (which exists in 5 dimensions) as it should. Thus the description by £ a is 
incomplete. 

However, several pages later he introduces the constraints 

£a.b,cA, e n c L d . e = 0. (9) 
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These constraints, (although not proven in Zhang's paper, since (^) is not granted) are the key for 
the effective interaction. If they are fulfilled, then one can conclude, that L is of the form 

L a ,b = c a n b - c b n a . (10) 

(This can be found in the following way: If rt\ = 1, and all other rii vanish, then one finds from 
(R) L a b = 0, if both a / 1 and 6^1. However L a \ — —L\_ a = c a for a ^ 1 can be chosen with 



arbitrary c a . Thus one obtains in this special case e q. (|lO|) . The relation (10) is form- invariant under 
SO(5)-rotations. Thus it holds for general rii.) Eq. (|10[) has two consequences: 

(i) One can now parametrize H a in terms of the c% and determine the minimum as one varies the 
coefficients c. One finds the minimum for 

2fin 5 2fi ni 

C 2 = C 3 = C4 = 0, Ci = — — , c 5 = — (11) 



with 



N= 1 J + 1 6 4 . (12) 

Xc Xir 



From this one obtains the effective interaction 

F e£F = -2/i — — -(n 2 +n 3 + n A ). (13) 

For Xc = X-n it agrees with V^ hang . If \c ^ Xx, then the effective potentials are different (actually in 
first order in Xc — Xt there is still agreement), but many conclusions Zhang has drawn will continue 
to hold. 

(ii) The constraints imply, that out of the representations (Ai,A2) the only allowed states have 
representations of the form (Ai,0). Let me explain this shortly. The irreducible representations of 
SO(5) are characterized by two (non negative) numbers (Ai, A2) with Ai > A2 > 0. The first number 
is the largest eigenvalue of one of the operators L, say L\ t i. Next consider the subspace of the states 
in this representation with this eigenvalue Ai for L1.2. Then A2 is the largest eigenvalue of one of the 
operators L, that commutes with the first one, which may be £3,4 in this subspace. Semiclassically, 
that is for large Ai and A2, one may replace the operators L a j, by their expectation values. Then the 
eigenvalues of the antisymmetric matrix L are ±iAi, ±iA2 and 0. If L has the form (|l0|), then 



Therefore the main problem left is to understand why the orthogonality constraints should hold. 
This will be done in the next section. 



3 Entropy as Source of the Orthogonality Constraint 

We claim that the orthogonality constraints are an effect of entropy. It will be shown that a system 
even without any interaction has an entropic parts of the form 

- (ta,b,cA,en c L dye ) 2 (15) 

which make it favorable for the system to obey the orthogonality constraints. Thus (0) is not a strict 
constraint, but if it is fulfilled, then the entropy assumes a maximum. 

3.1 Entropy of the Heisenberg Antiferromagnet 

First consider the Heisenberg antiferromagnet as an example, which Zhang has rightly mentioned in 
his paper. The Heisenberg antiferromagnet consists of two sublattices with magnetizations mi and 



3 



rri2. We expand the entropy in powers of m 2 and m|. Since the two sublattice-magnetizations describe 
the behaviour on different sublattices, the entropy is a sum of contributions on these sublattices 

S = S - Cl (m? + m 2 2 ) - c 2 ((m 2 ) 2 + (m 2 ) 2 ) - ... (16) 

Then denoting the homogeneous magnetization by mo and the staggered magnetization by m ?0 we 
may write mi = mo + m go and m 2 = mo — m 9o and obtain the entropy 

S = S - 2ci(mQ + m 2 qo ) - 2c 2 (mo + m 2 J 2 - 8c 2 (m • m 9o ) 2 - ... (17) 

Thus if there is no coupling in the interaction between mo and m go , the system prefers to have 
the homogeneous and the staggered magnetization orthogonal to each other due to the contribution 
— 8c 2 (mo • m 9o ) 2 in the entropy. Remember that the entropy enters into the free energy with a minus 
sign (F = E — TS) , and thus m • m qo — will yield the minimum of the free energy. Thus there is not 
a strict constraint m • m qo = for the Heisenberg antiferromagnet, but there is a term proportional 
to the square (m • m go ) 2 in the free energy, which favors the constraint to be obeyed. 



3.2 The entropy 

It will now be shown that the entropy of the more general SO(5)-invariant system contains terms of 
type (e a ,b,c,d,en c Ld.e) 2 ■ We will actually enlarge the system of operators L a ^ to those of an SO(8) (for 
details see in the appendix). To determine the entropy we start from the Hamiltonian 

Hq, = 2^ ^a : bL a ,b, Qb,a = —^a,b- (18) 
a, 6 

where the set of our L a j, also includes the n a — L a $. We do no longer use the normalization (|J). 
We note, that recently a U(4)-scheme including the operators L a ^ with a, b = 0..5 and Lq^ and 
their sub-groups has been considered in [^| (note that SU(4) is isomorphous to SO(6)). The f2 are 
introduced as Lagrange multipliers and will be adjusted to yield given expectation values of L a ^, and 
the entropy will be calculated up to fourth order in L. 

fi is an eight-dimensional antisymmetric real matrix. Since the matrix is antihermitean, its eigen- 
values are purely imaginary and occur pairwise. Thus the eigenvalue equation can be written 



E 



n a , b (x { b k) ± iyi k 3 ) = ±iJV( X W ± iy (*>) (19) 



with k — 1..4 and real vectors and y^ k \ The vectors and y( k > are orthogonal to each other. 
If they are normalized, then f2 may be represented 

n^ = E w( W-W). (20) 

k 

Next we perform a special orthogonal transformation, so that x^ and are oriented in appropriate 
directions, e.g. so that x^ x \ point in the 5- and 1- direction, x^ 2 \ y^ in the 2- and 3-direction, 
x^ and ?/ 3 ) in the 0- and 7-direction and x^ and j/ 4 -* in the 4- and 6- direction, resp. After this 
special orthogonal transformation H reads 

H tiaus = 2(o;( 1 )< 1 , 1 + a;( 2 )iV iliCT ,+c,W< i9il -c,W< 9i(r ,) 

= Y,^ il) +^ 2)s + ^ 3) 9( k )-^ i) 9(m(cl s Ck^^\). (21) 

k,s 

(with s = ±1, g(k) = ±1, g(k + qo) = —g(k)). Thus for momenta k and k + qo one has in total 2 4 
states constructed out of four single-particle states, which may be either occupied or unoccupied and 
which contribute the energies ± \t s ,g depending on whether the state is occupied or unoccupied 



u W+ u V>8 + uiWg-wWg8. (22) 
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Then we obtain for the partition function Z 

2 v ' 2 ^ 2 12 ' 2 

Summation over the four states yields 



InZ = ^ ln(exp(l/3e) + exp(-±/?e)) = ]T(ln2 + ^ - ^ + O^ 6 )), z = ^. (23) 



$> 2 = /3 2 ]>>« 2 (24) 

i 

E* 4 = ^(E wW4+6 E wW2w °' )2 - 24tj(1)w(2)w(3)w(4) )- ( 25 ) 

i i<j 

Similarly we expand the entropy S 

S = k B J2(^2-^ + ^+0(x e )) (26) 



and determine the expectation value of the quantity ZiM conjugate to (only the first component 
is given; the others are obtained by permutation), 

= J^iL = ^(D _ ^ ( w (l)3 + 3^(1)^(2)2 + J3)2 + w (4) 2) _ ^(2)^(3)^(4)) + q^S) (2?) 

and express w' 1 ' in terms of iW, 

/3a/ 1 ) = 2L« + ^(i^ 3 + 3L«(i( 2 ) 2 + L< 3 ) 2 + L< 4 ) 2 ) - 6L< 2 >L( 3 >L< 4 >) + 0(L 5 ). (28) 

In diagonal form one has = £5,1, etc. Then S for the subspace of the electrons with momenta k 
and k + qo reads 

5/fc B = 4 In 2 - 2 V i« 2 - - ( V L« 4 + 6 V L^ 2 L^ 2 - 24L^L^L^L^) + 0(L 6 ). (29) 

3 

% i i<j 

Thus we have expressed the entropy in terms of the eigenvalues ±iL^ of the 8 x 8-matrix (L a ^) of 
the expectation values L a _ b . We will express it now by the matrix-elements of L. First we have 



E i(i)2 = -^ 2 = -E^< (3°) 

E i(iH = = E L U + 2 E L lb L lc + 4 ££.,c£m£6,c£m. (31) 

z a<6 a.b<c a<b,a<Cc<Cd,b^c,b^d 

]>>« 4 + 2]>>« 2 L« 2 = 1 -(trL 2 ) 2 = ^Ll b + 2^L 2 , fc i 2 , c + 2 £ L 2 , b L 2 d , (32) 

i i<J a<fc a,b<c a<b,a<c<d,b^c,b^d 

L (l) i (2) i (3 )i (4 ) = pf(i)) (33) 



where pf (L) is the Pfaffian of L. In our case the indices a, 6 of the elements L a ^ are numbered from to 
7. Then the Pfaffian (which is defined only for antisymmetric matrices) is the sum Ylk=i ( — ) fe_1 -^o,fc-Pfcj 
where P fe is the Pfaffian of the matrix obtained from the matrix L by deleting the rows and columns 
with index (that is the first one) and index k. One continues recursively until the Pfaffian of a matrix 
with no entry is left, which is defined to equal 1. (If one starts with a matrix of odd dimension, then 
finally one arrives at the Pfaffian of the lxl matrix with entry 0, since it has to be antisymmetric. 
This Pfaffian is defined to equal 0. Therefore the Pfaffian of a matrix of odd dimensions vanishes.) 
We mention that the determinant of an antisymmetric matrix equals the square of its Pfaffian. For 
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2x2 and 4x4 matrices the Pfaffians read 



Pa,b,c,d ■= pf 



L„ 



L a ,bL c ,d ~ L a . c Lb.d + L a rfLf, , 



Then we have 



, ( L a , 6 

/ L a ,b L a .c L a ,d \ 

— L a ,b Lb^ c LbA 

—L ac —Lb.c Lc.d 
\ —L a 4 —LbA —L c .d J 

S/k B = 41n2 + trL 2 - i(|(trL 2 ) 2 - trL 4 - 24pf(L)). 



(Note that trL 2 is negative). Further algebraic manipulations yield 
-trL 4 = 4(trL 2 ) 2 + 4 £ p 2 ,^, 

a<b<c<d 

24pf(L) = 4 2J £a,b,c,d,eJ,g,hPa,b,c,dPe,f,g,h- 

a<b<c<d, e<f<g<h 

Putting all the contributions together we obtain 



(34) 
(35) 

(36) 

(37) 
(38) 



S/k L 



4 In 2 + trL 2 
4 



1 

12 



2 ,2 



(trL 2 ) 



^ (P0,a,6,c — ^ e a ,b,c,d,e,f,gPd,e,f,g) + 0(L 6 



(39) 



0<a<fc<c 



d<e<f<g 



This corresponds to the separation of (|2^) into 
S/fc B = 41n2-2^LW 2 -i(]TL« 2 ) 2 

_ !((i(DL( 2 ) - L^LW) 2 + (L«L( 3 ) - L' 2 'lW) 2 + (L«L< 4 ) - L^L^) 2 ) + 0(L 6 ). (40) 
3 

Thus the entropy consists of a completely rotational invariant contribution depending only on trL 2 
and a negative sum of squares of Po,o,&,c ± Pd,e,f,g- Note that the sum over d, e, /, g contains exactly 
one term ±p. Thus a maximum of the entropy is reached, when the arguments of all the squares 
vanish. Provided L a b vanishes if a or b equals 6 or 7 as assumed in the SO (5) theory, then out of 
the 35 squares 20 vanish identically, 10 have the form e a ,b,c,d,en c Ld, e , which when required to vanish 
are Zhang's orthogonality constraints, and 5 have the form £ a ,b,c,d,eLb,cLd,e- One can easily see from 
(|l0|), that if Zhang's orthogonality constraints are fulfilled, then also these latter quantities vanish. 
Thus the orthogonality constraints are not a strict requirements, but their fulfillment lowers the free 
energy. 



We mention, that the vanishing of the squares in the second line of (40) implies two types of 
solutions, either one L^ can be different from zero and the other L^ vanish, or all L^| are equal, 
and the product of the four L^' is positive. 

At half-filling (no doping) one has L\£ — 0. Then in the antiferromagnetic state, that is for non- 
vanishing components L2.0, L^^, and L^o, but otherwise vanishing components L, only one eigenvalue 
LW is different from zero. As soon as one has some doping L1.5 is different from 0. If the system is 
still antiferromagnetic, then more than one eigenvalue differs from zero. On the other hand, if the 
modes associated with n§ = Lq^ or nj = L-j^, are massive, as assumed, then it costs energy to have 
all four eigenvalues LW different from zero. Then it becomes preferable to have the superconductiong 
components n\ = L1.0 and n§ = L^^ different from zero and to have vanishing antiferromagnetism, 
since then again one has only one nonvanishing eigenvalue L". 
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4 Limits for the ground state 



The entropy argument given in the preceeding section can be applied in the vicinity of the critical 
temperature. At low temperatures the contribution of the entropy to the free energy decreases, since 
it enters with the factor T . Therefore we consider separately the situation at low temperatures. 



4.1 Antiferromagnet 

Let us start with a simple consideration for the antiferromagnet. Assume again two sublattices with 
magnetization mi and m 2 . Assume they are restricted by upper bounds |mi| < 1, |m 2 | < 1. Suppose 
now the system has a homogeneous magnetization mo. If now a staggered magnetization parallel to 
mo is added, then apparently |m go | < 1 — |mo|. If however the staggered magnetization is perpen- 
dicular to the homogeneous one, then one has the weaker restriction |m go | < y/l — m§. Thus the 
antiferromagnetic interaction can act more strongly, if the staggered magnetization is perpendicular 
to the homogeneous one. Thus again mo • m go — is fulfilled. 



4.2 Hartree-Fock-Bogoliubov ground state 

We have seen, that at low temperatures the bounds on the appropriate quantities (order-parameters) 
are important for the ordering of the ground-state. Therefore we will consider the bounds of the 
eigenvalues ±iL^ of the matrix L. If for a fixed matrix f2 in ([l8]) one takes the low temperature 
limit [3 — > oo, then normally a pure Hartree-Fock-Bogoliubov state instead of a mixed state remains. 
Depending on the sign of e s , g the occupation number n s , g assumes one of its extremal values and 1 
in the diagonal representation iTtrans 

n s , g = -- ^signe s , s . (41) 
From this we obtain the eigenvalues £,W 

L {1) = i ]>>,,-! (42) 

s,g 

L (2) = \T, sn s, 9 (43) 

s,g 

L {3) = \T.9n s , 9 (44) 

s,g 

L {4) = \j2 s 9 n ^- (45) 
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Putting now n s g — or 1 in all combinations one finds two types of solutions for L^ l \ In the first 
class one L^ 1 ' — ±1 and the other L^' = 0, whereas in the second class one has L« = ±\ for all % 
with the restriction for the signs, that the product of all lS l > is positive. One easily realizes that in 
all these cases the squares in the second line of eq. (^) vanish. As a consequence the orthogonality 
constraints are fulfilled for these states. The eigenvalues obtained for the Hartree-Fock-Bogoliubov 
ground-state are the extremal ones. For correlated states these eigenvalues can only be reduced. 
More precisely, the range of the set of possible eigenvalues ... lies in the convex volume 

bounded by the extremes given above. For an interaction quadratic in the operators L the energy 
assumes extremal eigenvalues for extremal eigenvalues L in the case of symmetry breaking. Without 
symmetry breaking the eigenvalues would vanish. We have argued before in favour of only one non- 
vanishing eigenvalue LW. Reduced to SO(5) this would be Ai, whereas the second largest (vanishing 
one) is A 2 of the representation (Ai, A 2 ). These correspond to the representations (Ai,0) actually found 
in the numerical calculations [|[ |). 
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5 Conclusion 



We have shown that the orthogonality constraints introduced by Zhang in his SO(5)-theory of highT c - 
superconductivity play an important role in the mechanism for the transition from antiferromagnetism 
to superconductivity as a function of doping. We have further shown that these constraints are not 
strictly fulfilled, but that their fulfillment yields a maximum of the entropy (for fixed J2^1b)- At low 
temperatures the entropy plays a weaker role. However, if the interaction drives the ground-state to 
the case of extremal eigenvalues which may be very well the case for an interaction bilinear in 
the operators L, then again the orthogonality constraints are fulfilled. 

We have expanded our scheme to a (mathematically natural) SO(8)-scheme. In this scheme two 
types of solutions for maximal entropy or extremal eigenvalues appear. Since the added degrees of 
freedom are probably massive, only those solutions, for which one eigenvalue is different from zero 
and the other ones vanish yield the minimum in the free energy. 

These considerations did not take the microscopic interaction seriously into account. Work in this 
direction has to be done. 

Acknowledgment I am indebted to Andreas Miclke and Jiirgen Stein for useful comments. 

A Operators in an SO(8) space 
A.l Bilinear operators 

All the operators in the SO(5) theory are of the form 

Kfn = ^E/W4 +g , s (7^kA t , (46) 

k,s,t 

</,7 = \ E /(*h.,t(4+,,.<*,t - l<oSs,t), (47) 

k,s,t 

= \ E f( k ) C k+q,s(v V l)s,tC-k,t- (48) 



Kin 



k,s,t 



Here the summation k runs over the Brillouin zone. The vector q can either be or q , where 2q is 
a reciprocal lattice vector. The staggered magnetization is described by the wave- vector qo as above. 
The function f(k) stands either for 1 or g(k) with g(k) = g(—k) = —g(k + qo) = ±1. Finally 7 is a 
hermitean two by two matrix. It may be either the unit matrix or one of the Pauli matrices. 
First we consider the symmetry of N + . If we exchange the two c '''-operators, then we obtain 

N qJ,l = N q,f(.+ q ),<jy 7 T a y (49) 

One has 

/(•+«>-«^/-»-^-{_^;:£^;i,. « 

and 

o v ~< T o v = s-v7 with 8~ = < . , T ^° r i 7 , . (51) 
' ' [ — 1 for Pauh matrices 

From this we conclude that only the six operators Nq 11: Nq lt N^ q11 , N+ g (Ja with s q jSj — 1 are 
different from zero. The same holds for N~. 
The hermitean adjoint operators are 

( N qjJ = S «-fKf^> ( 52 ) 
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!/ t 




71-t = 
" a 


N+ a 


T ~ l f t 
i>5,l+a = -y^ot _ 




7T Q = 




Ll+a,l+/3 = e a,/3,7 


"50,7; 


So,-y : 


- iV 0,l,erT ) 


ni = Li |0 = i(A f 4 


-A), 


At = 




ns = £5,0 = ^(A f - 


-A), 


A = 


~ iN 0,g,l> 



We obtain for the commutators 

[^T/,7'%"/',y] = o. (55) 

\. N qJ,l^ N q',f',l'] = 2^+«'' / (- + «') / '' Ty ' _iV 9+«'.//'(-+«).7'7)- ( 59 ) 

A. 2 Zhang's Operators 

The operators introduced by Zhang, which obey the commutator relations (|^) of an SO-group were 

£5,1= Q =JVj,i.i. ( 6 °) 

(61) 

(62) 
(63) 
(64) 

(65) 

n-l+a = ^1+q,0 = Sqa,a = N q0tltaa . (66) 

Here we have added n a = L a ,o to the group, since they obey the same commutator relations. 
A. 3 Extension to operators obeying an SO(8) group 

Obviously one can expand the range of operators L a j, by including the other operators N introduced 
above so that any pairs of particles, of holes or particle-hole pairs with total momentum or qq in 
the singlet and triplet channel appear. This allows us to introduce components L a jj with a, b = 6 or 
7. This expansion does not imply, that the symmetry group for the phase transition will be enlarged. 
Actually one expects that the new components of the 'order-parameter' stay massive at the transition 
and will not contribute directly to the symmetry-breaking. 

It is wellknown ||, [l0| that there are additional operators adding or removing two electrons, 
and r], which commute with the Hubbard Hamiltonian, which allows us by means of the commutator 
relations (|^) to introduce 

£1,6 = ^-^), ^ =<,!,! ( 67 ) 

£5,6 = ^+77), V = N- gAA , (68) 

Li+as = -NS^, (69) 
n 6 = L 6)0 = iN° ^ v (70) 

Finally there are 7 more operators left, which fulfill the appropriate commutator relations with 

Li,t = \tf+fi), i? = N+ ltl , (71) 
L5,r = l(tf-V), fj = No, hl , (72) 
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-^6,7 

n-7 = Lr,o 



iN° „„„ 
9o,g,c-° 



= -N, 



90,1,1' 

-AT 



(73) 
(74) 
(75) 
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